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Abstract 

An analysis of experimental data from the RFX (Reversed Field experiment) 
reversed field pinch machine [L. Fellin, P. Kusstatscher and G. Rostagni, Fu- 
sion Eng. Des. 25, 315 (1995)] is carried out to investigate the possible 
existence of deterministic chaos in the edge plasma region. The mathematical 
tools used include Lyapunov exponents, Kaplan- Yorke dimension, minimum 
embedding dimension estimates and nonlinear forecasting. The whole analysis 
agrees in ruling out the possibility of low-dimensional chaos: The dimension 
of the underlying dynamical system is estimated to be > 10. From a critical 
re-reading of the literature it emerges that the findings of this work are likely 

to be common to all reversed field pinches. 
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I. INTRODUCTION 



Fusion machines are complex systems possessing a huge number of de- 
grees of freedom where phenomena involve many scales of length and time. 
Such devices are therefore natural candidates for checking of nonlinear theo- 
ries. In the past years, as soon as the experimental instrumentation made it 
possible, some effort has been exerted to determine if the temporal behaviour 
of measured quantities in fusion plasmas, displaying a fluctuating behaviour, 
could be modelled using some low-dimensional deterministic theory or were 
the manifestation of truly random processes . The results of this search 
for deterministic chaos are contrasted: while the results of numerical simu- 
lations show the existence of low-dimensional chaos, experimental works 
find positive (2|-|5| as well as negative results j?]|§. 

In this paper we present a study on this subject done on experimental 
data taken from the RFX experiment: a large (R = 2m, a = 0.457m) toroidal 
device for the magnetic confinement of plasmas in reversed field pinch (RFP) 
configuration Q built in Padova (Italy) and designed to reach a plasma cur- 
rent of 2 MA. The data refer to measurements made in the outer region of the 
plasma, inside the last closed magnetic surface. A number of mathematical 
tools has been applied to these data to detect traces of deterministic chaos. 
Furthermore, a critical discussion of the existing literature on the subject 
is done, in the light of the results found, and some general conclusions are 
drawn. 

In the next section we illustrate the mathematical apparatus used, then a 
brief description of the diagnostics used to collect the data is given. Section 
IV presents the results of the analysis and finally some comments are given. 
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II. NUMERICAL TECHNIQUES 



Several analysis techniques have been developed to identify the presence 
of a low dimensional attractor in a dynamical system from its time series. In 
this section we describe the mathematical tools used in this work. Some of 
them are very recent and only recently have been applied to real experimental 
data. Therefore there is a further interest in verifying how they perform in 
this situation. All the techniques presented herein apply to a scalar time 
series of recorded data S = (x±, X2, ■ ■ ■ , xn). From the scalar series S one 
may construct g-dimensional vectors in delay-coordinates, 

= [ x ii X i+T, ■ ■ ■ i x i+(q-l)r] 

The (/-dimensional space of delay vectors plays a central role in this kind of 
analysis. The integer r is known as the time lag. It can be larger than 1, i.e. 
data may be sampled at a frequency lower than the experimental one. This 
is done since choosing r smaller than the autocorrelation time t c of the data 
would introduce spurious correlations. On the other hand r should neither be 
too large. The standard choice is to put r ~ 2-\r3t c . In our units, depending 
upon the data, t c = 1 -r 2. To test the reliability of our results, r was varied 
in the test from 1 to 7 times t c . Conclusions were unaffected by the choice. 



A. Estimate of the minimum embedding dimension 

In experiments, the true dimension of the phase space of the system under 
study (the embedding dimension m) is usually unknown. A reasonable guess 
for it is however essential for any analysis. It has been proven by Takens (10] 



that for a chaotic dynamical system the time recording of a single variable is 
sufficient to reconstruct the relevant dynamics-and in particular the dimen- 
sion d of the attractor-provided that m > 2d + 1 (in actual cases it is known 



that this condition may be slightly relaxed). 

In this work we have implemented the method proposed by Cao [11| to 
estimate the minimum embedding dimension from a scalar time series. 
This method uses the time-delay vectors (with r set to one for simplicity) 

Yi{m) = [xi,x i+1 ,x i+ 2, . . • , x i+(m _ 1} ] , i = 1, 2, . . . , N - (m - 1) 

(where m is the guessed embedding dimension) to build the function 

f . v \\yi(m+ 1) -y„(i, m )(m + 1)|| . 

a(z, m) = n — - — '— — , i = l,2,... ,N -m, 

||y;(m) - y„( iim )(m)|| 

where || • || is a norm and n(i,m) is an integer such that y n (i^ m ){ m ) ls the 
nearest neighbour of yj(m) in the m-dimensional space, according to the || • || 
norm. The actual functional form of 1 1 • 1 1 does not appear to be of importance. 
We have used the maximum norm: 

||Y — Z| | = max \ Yi — Zi\ , i = 1, . . . , m 

If m is the true embedding dimension, any two points which are close together 
in the m-dimensional reconstructed space, will stay close also in the m + 1 
space. Points which satisfy this condition are called true neighbours, otherwise 
false neighbours fl2|| . Starting from a low value for m and approaching the 
correct value (hereafter referred to as v) the number of false neighbours should 
decrease to zero, or equivalently a(i, m) should reach a constant value. Cao 
p| suggests to use averages of this quantity: 

N—m 

E ( m ) = T7 E a (h m ) 

N — m f— f 

i=i 

and 

El(m) = E(m + l)/E{m) , 

which allow to obtain results independent upon the sample data chosen. 
£71 (m) should stop changing when m becomes greater than some value v, 



if the time series has a finite dimensional attractor. In the case of random 
data El will never saturate but when dealing with real data it is difficult to 
distinguish if it has attained a constant value or is slowly increasing, therefore 
in Ref. [11| it is recommended to also compute the function 



E2(m) = E* (m + 1) / E* (m) , 

-, N-m 

E*{m) = -^— J2 \x i+ m-x n{itm)+m \ , 
i=i 

where the meaning of n(i, m) is the same as above. For random data E2 will 
stay close to 1: the x's are now independent random variables, and therefore 
their average distance will be the same regardless of the space dimension 
m. For deterministic data, there will be a certain correlation between them 
which makes E2 a function of m. As an illustration of the method, we plot 
El, E2 in Figure 1 for two sets of data: one time series is computed using the 
Mackey-Glass equation 

dx(t) , . 0.2x(t - A) 

-O.lx(t) + 



dt ' w l + x(t- A) 10 

with A = 30, which is known to describe the dynamics of a chaotic system 
with an attractor dimension of about 3.6 [13]. The other time series is gen- 



erated using a random number generator. In the random data series E2 is 
always very close to one and El slowly converges to the same value. The plot 
of the deterministic map shows instead that both El and E2 reach the same 
value at D ~ 6; furthermore the behavior of El is not that of an asymptotic 
convergence to 1, but is more akin to the reaching of a threshold value. 



B. Lyapunov exponents 



The Lyapunov exponents, measuring the average divergence or conver- 
gence of orbits in phase space, are among the most frequently used quantities 



to ascertain the presence of chaos. Given a map x(t) = f xo (^)) f : R m ~~ > R"\ 
the Lyapunov exponents are defined as 

\W = limiln|^#|= lim±||J t (fc) ||, k = l,...,m , 
(k) 

where ||J t || is the kth eigenvalue of the m-dimensional Jacobian. A nec- 
essary condition for a system to be chaotic is to have at least one positive 
exponent. 

Several well established techniques exist to compute the Lyapunov expo- 
nents for a system whose dynamical evolution is analytically known. Extract- 
ing them from an experimentally determined data set is much more difficult 
due to the limited length of the sample and to the presence of noise. Ex- 
isting algorithms usually fit the experimental points to an analytical map 
g : R m — > R such that 

g ([xj, . . . , Xi-\- m — l]) = %i+m 

or, which is the same, to the map 

g([xi, . . . ,x i+m -i\) = [x i+ i, . . . ,x i+m ] , g:R m ^~R m . 
Under quite general conditions, the largest v Lyapunov exponents of f and 



g are the same |1JJ15|. Therefore one estimates the A's using standard tech- 
niques on g,g. In this work we have used two codes for estimating Lyapunov 
exponents. The former is the code developed by Watts HR], based upon the 



method by Briggs [14,nj, and already used on data from a magnetically con- 
fined plasma. In this method the time series is embedded in a delay space of 
given dimension, a number of nearest neighbours is found, and their trajec- 
tory is fit to an analytical function (usually a polynomial: in our runs we have 
used a polynomial of order 2). Then the Jacobian (10) may be obtained by 



analytical differentiation. The second code has been developed by one of us 
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17] based upon a method by Gencay and Davis Dechert jjOj]. Here, a single 
global fit is done by using logistic maps, i.e. 

L L n 

i 1 + exp(-bi - w ; ■ Xi) 
where L is the number of functions v and each w is an m-dimensional array. 
The Vs, A' s an d w z' s are fitting parameters. Logistic maps have some inter- 
esting features: they may fit arbitrarily well any analytical function as well as 



its derivative. In Ref. [15| it is stated that this choice of functions has some 
advantages over the local polynomial approach in terms of stability of results 
in presence of noise and of fewer needed data points. On the other hand it 
requires a nonlinear fitting which is computationally more demanding. 

C. The correlation and Kaplan Yorke dimensions 

The correlation dimension D c gives another estimate of the embedding 
dimension of the system. D c is defined as 

1 N 

Cd{r) = li^ooN! £ B(r-\Xi-*i\) , 

v£j=l 

where H is the Heaviside function and the X's are defined in Eq.(jl]). For 
m > v, lim r ^o C(r) ~ r v . 

The correlation dimension has been and currently is a favourite tool to diag- 
nose the presence of chaos in fusion plasmas |||§. 

Instead of D c we have computed, using Watts's code the Kaplan- Yorke 
(or Lyapunov) dimension 

L>KY — J H : 

-Aj+i 

where j is the largest integer such that J2i=i > 0, with the A's ordered 
as Ai > A2 > . . .. There exist some conjectures [^] according to which the 
uguagliance Dky ~ D c holds 



Dky is much easier to compute than D c , since it is only necessary to know the 
Lyapunov exponents, previously computed. This strength is at the same time 
a weakness when there may exist uncertainties about the correct value of the 
exponents. This is our case, however our analysis does not rely just on this 
single parameter and as we shall see, all results corroborate the indications 
from Dky- 



D. Nonlinear forecasting 

Given a time series of finite length representative of a chaotic dynamical 
system governed by the map f : R m — > R m , an inverse problem consists in 
finding a smooth map g : R m — > R m , or its projection g : R m — ► R such 
that g be an accurate approximation of f. In the case g is obtained, it may 
be used to predict accurately further data points. Otherwise, if the system 
is not governed by a finite dimensional map, or if the guessed m is too low, 
the forecasting will be unreliable after few predictions. A review about the 



subject may be found in Ref. [Jl 3 1 . 



We have used as fitting function the logistic maps (Eq. |13J): the original 
time series S of length N has been divided into two parts of lengths N\ and 



.ZV2 = N — N\. The first JVi data have been used to fit Eq. (13) and the 
predictive power has been tested on the remaining N% points. The measure 
of goodness is the predictive error: 



2 1 ( x n di^-n—m)) 

a =W„ 2- ~2 



-^2 „ , i <*i 

n=Nj+l - 1 

where a^. is the variance of the time series (in this work all data have been 
normalized so to have a^. = 1). The increase of the trial embedding dimension- 
provided that enough fitting parameters are allowed-will give a slight decrease 
of a 2 until m ~ v, when a sudden decrease to much smaller values is expected. 
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To provide the reader with an example, in Figure 2 we have plotted a 2 versus 
m for the Mackey-Glass map where is clearly visible the decrease of more 
than two orders of magnitude of a when the dimension is > 4. For comparison 
we have estimated the predictive power of the method against a time series 
generated from a normal distribution. As expected the possibility of any 
forecast is null, with a 1 always greater than the variance of the original data. 

It is worth mentioning that estimating the Lyapunov exponents and corre- 
lation dimensions has been attempted in all of this kind of studies concerning 
fusion plasmas. The forecasting of the data, however, has been applied, to 



our knowledge, only to data from the Madison Simmetric Torus (MST) [19] 
reversed field pinch device Qj. The estimate from the time series itself of the 
embedding dimension, finally, is applied in this work for the first time. 

Some words must be spent about the confidence which may be assigned to 
the algorithms. Two crucial topics affecting their performances are: the num- 
ber of data available and their quality, i.e. how much they are polluted with 
noise. As far as the first point is concerned, usually the more data one can 
elaborate the larger may be the dimension of the system which may be cor- 
rectly estimated. We had available records of some thousands values, about 
as many as used by MST group. The authors of all algorithms used here claim 
them to be able to reconstruct the correct dynamics of a low-dimensional sys- 
tem (i.e. four or five-dimensional) using few hundreds data, with the possible 
exception of the Watts' code (see [||,|7]]). In conclusion, and by comparison 
with the MST group's estimates, we may assert that dimensions up to or just 
below 10 may be correctly detected by our techniques. 
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III. RFX EXPERIMENTAL DATA 



Signals coming from three different diagnostic techniques have been anal- 
ysed in search of low dimensional chaos features: 

(a) Floating potential (Vf): This is the potential of an electrically insu- 
lated conducting probe immersed in the edge plasma. It is known to be related 
to the local plasma potential V p and to the local electron temperature T e (in 
eV) through the relationship Vf = V p — aT e where a is a constant which 
for the RFX edge plasma is approximately equal to 2.5. The probe was a 
graphite pin housed in a boron nitride structure which protected it from the 
unidirectional superthermal electron flow commonly observed in RFP edge 
plasmas. Data were sampled at 1 MHz. The measurements were collected 
during the experimental campaign described in [^0) . 

(b) Time derivative of the radial magnetic field (dB r /dt): It was measured 
with a pick-up coil housed in the same boron nitride structure as the Vf 
measuring pin. This is a local measurement, like the previous one. Data were 
sampled at 1 MHz. 

(c) Density fluctuations at the edge: Collected by a reflectometer. The 
RFX reflectometer has been especially designed to deal with high-frequency 
fluctuations. It is an homodyne reflectometer |^] operating at a maximum 
sweep rate of 4 GHz///s in the range 34-38 GHz. The collected signal is of 
the form 

s(t) = A(t) cos(Ac^)) 

with A amplitude and A(f> phase difference of the reflected radiation at a fixed 
microwave frequency. In this work we have analyzed the temporal behaviour 
of the amplitude. 

Figure 3 displays the plasma current waveform and typical waveforms for 
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the three signals considered. The data were collected in discharges having a 
plasma current ranging between 350 and 600 kA. 

IV. RESULTS 

Figure 4 displays the quantities El, E2 (Eqns. || [?]) for the three analyzed 
signals. The data are taken for three different shots (shot 8422 for magnetic 
fluctuations, shot 7999 for potential measures, and shot 7852 for reflectometer 
data). A number of data points ranging up to about 10000 has been used. The 
resemblance with random data (Fig. 1) is impressive, which clearly suggests 
the existence of a very high dimensional phase space. 

A further confirmation is obtained by plotting the Kaplan- Yorke dimen- 
sion (Fig. 5): no sign of saturation is obtained and over the whole explored 
range Ojfy(m) ~ m. 

In Figure 6 we plot the predictive error (Eq. [[§). The value of a is always 
close to one, meaning that no real accurate prediction may be done. Some 
very small differences may be seen among the three signals, even if they may 
well be just subjective impressions: the magnetic fluctuations data show the 
worst predictive power, suggesting perhaps that some different mechanism 
from the other two diagnostics is at work. This could be the case if mag- 
netic fluctuations are mostly due to magnetohydro dynamics tearing modes 
whereas potential and density fluctuations are mainly caused by some elec- 
trostatic instabilities localized in the edge region. For the numerical aspect, 
the calculation of this quantity turned out to be rather reliable, with limited 
variations between runs. 

Finally, in Figure 7 we plot the largest Lyapunov exponents. Since the 
system does not appear to be dominated by low-dimensional chaos, the use- 
fulness of these coefficients is rather limited, however it may be interesting 
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to compare the two approaches. The estimation of the A's from real data is 
known to be a difficult task; we did not obtain stable results varying the pa- 
rameters: indeterminacies of 50% are quite likely, so the values shown are to 
be considered as representative of the general trend. We found that Watts's 
code is rather sensitive to the choice of the input parameters such as the order 
of the polynomial and the number of neighbours to be used in the fitting. Our 
code is more stable from this point of view since it needs very few input data 
(essentially, the embedding dimension m and the number of logistic map L). 
However, it is difficult to perform the fit over a large number of data points 
(600 is the maximum used), so the results suffer of the scarce statistics. 

In Watts's work it is emphasized that, in order for the previous analysis 
to hold, the physical system must be in a stationary state. The breaking 
down of this assumption may translate to an overestimate of the attractor 
dimension or even make it non measurable. Watts quotes two possible causes 
for the lack of non-stationarity: (a) The plasma may not reach a state of 
equilibrium; (b) Even if equilibrium is reached, random perturbations (Watts 
cites as an example the influx of impurity ions from the walls) may destroy 
it. In RFX discharges true flat-top periods lasting some tens of milliseconds 
are reached. In our study we have considered both discharges where this flat- 
top period was reached, and others where instead plasma parameters were 
slowly changing. Conclusions are unaffected by the discharge chosen. We 
are therefore confident that any lack of stationarity due to this causes may 
just very slightly modify the results. Point (b) is by nature uncontrollable: 
however, since the same results are obtained from different discharges it may 
be hoped that these random events have not affected the final results. 

All our results agree in clearly pointing out that no low-dimensional chaos 
appears at the RFX edge. This is well consistent with findings of Watts et 
al. Hf?! for MST. When limiting to RFP's the only other research within this 
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field is that done in HBTX1A ]22j by Gee and Taylor ]2j, where traces of 
low-dimensional chaos were found in magnetic-field oscillations. This sharp 
discrepancy puzzled us, therefore we resorted to check conclusions of Gee and 
Taylor against their own data. In Ref. Q the attractor dimension is estimated 
using the correlation dimension technique (Eq. |14|); their figure 1 shows the 
behaviour of log C (r) versus log r . Gee and Taylor claim that the slope of these 
lines saturate in correspondence of a dimension ~ 7. We could not check this 
since in their figure 1 only dimensions from 1 to 6 are plotted. However, from 
the data available, we could not find any clue of such a saturation: the slope 
of the straight line for m = 6, as estimated by visual inspection of the plot, 
appears close to 8. Even allowing for the large error induced by our gross 
way of estimating, the true result cannot be much smaller than 6, which is a 



necessary condition to speak about a saturation of the slope 23]. 

Up to this point, our work has focussed almost entirely on RFP plasmas. 
Some interesting considerations may be drawn by comparison with those toka- 
mak plasmas where the fingerprints of deterministic chaos have been found. 
We refer to the works §{§. From Ref. §, using data of TOSCA @] and the 
Joint European Torus (JET) |25| tokamaks, it appears that a small value of 
the dimension of the dynamical system is more likely to be found if the level 
of turbulence (measured for example by 8B/B) is small, which is not what 
happens in RFP's. Barkley's data [|| are a bit difficult to interpret since their 
analysis are done on filtered data, i.e. by selecting the wave number com- 
ponents k. Their finding is that the dimension increases with k. However, 
Barkley's data are chord averaged, so the central plasma plays a dominant 
role. In our work, and in the others studied, only edge quantities have been 
considered. It is quite possible that different mechanisms be at work in the 
two zones. 
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V. SUMMARY AND CONCLUSIONS 



In this work four statistical tools have been applied to the signals of some 
plasma turbulence measurements (magnetic field-, electrostatic potential-, 
and density fluctuations) of the RFX experiment, addressing the question of 
the existence of low dimensional deterministic chaos in them. The methods 
of analysis are well established (correlation dimension) as well as more recent 
(predicting errors and minimum embedding dimension estimates). 
All the conclusions are strikingly in accordance in ruling out that the dy- 
namics of the edge plasma in RFP's may be significantly affected by any 



low-dimensional process [26|. An estimate of a lower bound for the dimen- 
sion of the system may be given by Figure 6, which shows that this value 
must be greater than 10. This conclusion is further strenghtened by a critical 
re-examination of previous results Q conflicting with ours, which has shed 
some doubt about their validity. 

ACKNOWLEDGEMENTS 

The authors are very grateful to Dr. C.A. Watts for allowing them to use 
his code, and to Dr. R. Cavazzana for providing the reflectometry data. 



14 



REFERENCES 

[1] M. Persson and H. Nordman, Phys. Rev. Lett. 67, 3396 (1991). 

[2] S.J. Gee and J.B. Taylor, in Proceedings of the 12th European Conference on Controlled 
Fusion and Plasma Physics, vol. II (European Physical Society, Geneva, 1985), p. 446. 

[3] A. Cote, P. Haynes, A. Howling, A.W. Morris, and D.C. Robinson, in Proceedings of the 12th 
European Conference on Controlled Fusion and Plasma Physics, vol. II (European Physical 
Society, Geneva, 1985), p. 450. 

[4] H.J. Barkley, J. Andreoletti, F. Gervais, J. Olivain, A. Quemenur, and A. True, Plasma Phys. 
Contr. Fus. 30, 217 (1988). 

[5] G. Strohlein and A. Piel, Phys. Fluids B 1, 1168 (1989). 

[6] C.A. Watts, Chaos and Simple Determinism in Reversed Field Pinch Plasmas, Ph.D. thesis, 
University of Wisconsin (1993). 

[7] C.A. Watts, D.E. Newman, and J.C. Sprott, Phys. Rev. E 49, 2291 (1994). 

[8] M.L. Sawley, W. Simm, and A. Pochelon, Phys. Fluids 30, 129 (1987). 

[9] L. Fellin, P. Kusstatscher and G. Rostagni, Fusion Eng. Des. 25, 315 (1995). 

[10] F. Takens, in Dynamical Systems of Turbulence, edited by D.A. Rand and L.S. Young, Lecture 
Notes in Mathematics, Vol. 898 (Springer, Berlin, 1981), p. 366. 

[11] L. Cao, Physica D 110, 43 (1997). 

[12] M. Kennel, R. Brown, and H. Abarbanel, Phys. Rev. A 45, 3403 (1992). 

[13] M. Casdagli, Physica D 35, 335 (1989). 

[14] R. Brown, P. Bryant, H. Abarbanel, Phys. Rev. A 43, 2787 (1991). 

[15] R. Gencay and W. Davis Dechert, Physica D 59, 142 (1992). 



15 



[16] K. Briggs, Phys. Lett. A 151, 27 (1990). 

[17] F. Sattin, Comput. Phys. Commun. 107, 253 (1997). 

[18] D.A. Russell, J.D. Hansen, and E. Ott, Phys. Rev. Lett. 45, 1175 (1980). 

[19] R.N. Dexter, D.W. Kerst, T.W. Lovell, S.C. Prager, and J.C. Sprott, Fus. Techn. 19, 131 
(1991). 

[20] V. Antoni, R. Cavazzana, D. Desideri, E. Martines, G. Serianni, L. Tramontin, Phys. Rev. 
Lett. 80, 4185 (1998). 

[21] R. Cavazzana, F. Chino, M. Moresco, A. Sardella, and E. Spada, in Proceedings of the 24th 
European Conference on Controlled Fusion and Plasma Physics, vol. I (European Physical 
Society, Geneva, 1997), p. 361. 

[22] H.A.B. Bodin, C.A. Bunting, P.G. Carolan, L. Giudicotti, C.W. Gowers, Y. Hirano, I.H. 
Hutchinson, P.A. Jones, C. Lamb, M. Malacarne, A. A. Newton, V.A. Piotrowicz, T. Shimada, 
M.R.C. Watts, in Proceedings of the 9th Conference on Plasma Physics and Controlled Nuclear 
Fusion Research, vol. 1 (International Atomic Energy Agency, 1983), p. 641. 

[23] Notice that also Watts (see Ref. ||) did a critical re-reading of most of the papers dealing 
with chaos in fusion devices, reaching the conclusion that the results shown in many of them 
suffered of an unappropriate elaboration of the data. Gee and Taylor's results, however, were 
not explicitly questioned by Watts. 

[24] K. McGuire, D.C. Robinson, A.J. Wootton, in Proceedings of the 7th Conference on Plasma 
Physics and Controlled Nuclear Fusion, vol. 1 (International Atomic Energy Agency, 1979), 
p. 335. 

[25] E. Bertolini, Fusion Eng. Des. 30, 53 (1995). 

[26] This is not exactly the same as asserting that no low-dimensional processes exist. Actually, 
the signals may be composed by a contribution coming from these processes together with 



16 



another due to turbulence. What we can say is that the latter is overwhelming. 



17 



FIGURES 



FIG. 1. Estimate of the minimum embedding dimension using Cao's method (section |II A ). 



FIG. 2. Diamonds, predictive error for the Mackey-Glass equation (Eq. pi) versus embedding 



dimension as estimated by the method of section II D; Squares, the same for a series randomly 
generated from a normal distribution. 

FIG. 3. Examples of the signals used. Smooth curves are the plasma current (I p ), the wildly 
fluctuating ones are the signals. Note that we are referring to three different discharges. 



FIG. 4. El, El (see section II A| ) for the three signals. Here and in all the following plots V 



stands for potential fluctuations; B, magnetic field; n e , density fluctuations. Diamonds, El; Stars, 
E2. 

FIG. 5. Kaplan- Yorke dimension for the three signals. 

FIG. 6. The largest Lyapunov exponent versus embedding dimension. Stars are results from 
Watts' code, diamonds from our code. 

FIG. 7. Predictive error versus embedding dimension. 
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